Passive scalar transport in peripheral regions of random flows 
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We investigate statistical properties of the passive scalar mixing in random (turbulent) flows as- 
suming weakness of its diffusion. Then at advanced stages of the passive scalar decay its unmixed 
residue is concentrated primarily in a narrow diffusive layer near the wall and its transport to bulk 
goes through the peripheral region (laminar sublayer of the flow). We conducted Lagrangian numer- 
ical simulations of the process for different space dimensions d and revealed structures responsible 
for the transport that are passive scalar tongues pulled from the diffusive boundary layer to bulk. 
We investigated statistical properties of the passive scalar and of the passive scalar integrated along 
the wall. Moments of both objects demonstrate scaling behavior outside the diffusive boundary 
layer. We propose an analytical scheme for the passive scalar statistics, explaining the features 
observed numerically. 
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I. INTRODUCTION 

Stochastic dynamics of such scalar fields, as temper- 
ature or concentration of pollutants, in random (turbu- 
lent) flows is of great importance in different physical 
contexts, from cosmology to micro-fluidics. If the back 
reaction of the field to the flow is negligible then the field 
is called passive scalar. We consider the passive scalar in 
random flows, where the flow velocity varies randomly in 
time. Theoretical examination of dynamical and statisti- 
cal properties of the passive scalar in random flows 
back to classical works of Obukhov and Corrsin 
where a phenomenology for the passive scalar statistics in 
turbulent flows was developed in spirit of the Kolmogorov 
scheme 0. Modern understanding of the passive scalar 
statistics in turbulent flows is reflected in Refs. 0-0], 
see also the books [sl-Hoj. The mixing problem for the 
passive scalar is investigated for chaotic flows as well, see 
the book of Ottino 



11 



An interesting example of a ran- 
dom flow is the so-called elastic turbulence, discovered by 
Groisman and Steinberg [l2[ in polymer solutions. Ob- 
servations of the passive scalar statistics in the elastic 
turbulence were reported in Ref. (l3j . 

In 90th, a series of theoretical works devoted to the 
passive scalar statistics was done in the framework of the 
so-called Kraichnan model where the turbulent flow is 
assumed to be short correlated in time and possessing a 
definite scaling. The works enabled one to reveal general 
features of the passive scalar statistics in turbulent flows 
including the so-called anomalous sca ling and intermit- 
tency fl#[l6| . see also the surveys fl7l Il8l| However, the 
approach implies spacial homogeneity of the flow statis- 
tics and is not, consequently, directly applicable to re- 
gions near walls. 

In this paper we investigate the passive scalar statis- 
tics in peripheral regions of a vessel where the developed 
(high-Reynolds) turbulence is excited. Speaking about 



the peripheral regions of the turbulent flows, we imply a 
laminar (viscous) sublayer formed near walls where the 
velocity field can be regarded as smooth, it varies on 
distances of the order of the thickness of the sublayer. 
However, the velocity remains a random function of time 
there. Some laminar boundary layer is characteristic also 
of the elastic turbulence [lj| . The passive scalar statistics 
in the peripheral region is determined by a complicated 
interplay of its diffusion and random advection in highly 
anisotropic situation caused by presence of the walls. 

We are interested in advanced stages of the passive 
scalar decay assuming that the Peclet or the Schmidt 
number is large. In the case the unmixed fraction of 
the passive scalar is located mainly in a narrow diffusive 
layer near the wall, thinner than the thickness of the pe- 
ripheral region [2fJ. Then the passive scalar transport 
to bulk goes through the peripheral region outside the 
diffusive layer. Just this peripheral region plays a cru- 
cial role in formation of statistical characteristics of the 
passive scalar transport. The same reasoning can be ap- 
plied to a stationary case related, say, to a permanent 
heat flow going from the walls through the periphery re- 
gion to bulk. Moreover, fast chemical reactions can be 
analyzed inside the same scheme, see Ref. [2l|. A the- 
oretical approach to the problem was developed in Ref. 
[22| . principal predictions of the theory were confirmed 
by experiment, see Ref. [23j . 

To check the theoretical predictions in detail, we con- 
ducted extensive numerical simulations of the problem 
based on Lagrangian dynamics of particles representing 
the passive scalar. Aiming to establish main qualitative 
properties of the passive scalar transport in the periph- 
eral regions, we focused on the 2d (two-dimensional) case. 
However, a big advantage of the Lagrangian scheme is an 
ability to extend the approach without essential troubles 
to higher dimensions. We performed simulations for the 
space dimensionalities d — 3 4- 5 to establish universal- 
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ity of the passive scalar statistical behavior in 2d and to 
reveal features characteristic of higher dimensions. We 
used a scheme with permanent injection of particles near 
the wall that produces statistically homogeneous in time 
statistics of the passive scalar. However, our conclusions 
are correct for the decaying well because of adia- 

baticity: events responsible for the passive scalar trans- 
port to bulk are much faster than the average passive 
scalar decay. 

The obtained numerical data can be used to compute 
averages characterizing the passive scalar statistics. First 
of all, we found moments of the passive scalar at different 
separations from the wall. The data show an existence of 
a well pronounced diffusive layer where the passive scalar 
is mainly concentrated, in accordance with the theoret- 
ical expectations formulated in Ref. (22J. Outside the 
diffusive layer the passive scalar moments demonstrate 
scaling behavior with exponents deviating from ones pro- 
posed in Ref. [Hf where diffusion was assumed to be neg- 
ligible outside the diffusive layer. We checked that the 
deviations are related to diffusion, indeed. The situation 
resembles the passive scalar statistics in the Batchelor 
velocity field on scales larger than the pumping length 
where diffusion appears to be relevant [24J, that corrects 
the diffusionless behavior examined in Ref. (25|. Next, 
we introduced the passive scalar integrated along a sur- 
face parallel to the wall. The diffusion along the wall 
drops from the equation for the object. We demonstrated 
numerically that outside the diffusive layer moments of 
the integral passive scalar have well pronounced scaling 
behavior, as a function of separation from the wall. We 
found the corresponding scaling exponents for moments 
with degrees n = 1-^6 in space dimensionalities d = 2-^5. 
The moments exhibit an anomalous scaling signalling 
about strong intermittency of the passive scalar statis- 
tics. 

The simulations enabled us to reveal objects lying be- 
hind the intermittency. The objects are tongues of the 
passive scalar pulled from the diffusive layer towards 
bulk. The tongue cross section diminishes as the sep- 
aration from the wall increases (due to increasing the 
velocity component perpendicular to the wall). That ex- 
plains why diffusion can play an essential role even in 
the region outside the diffusive layer. The subsequent 
tongue evolution including tongue folds produces long- 
living structures possessing complex shape. Sometimes 
the tongues are pulled so strongly that push irreversibly 
a passive scalar portion away from the wall. Just this 
mechanism is responsible for the passive scalar transport 
to bulk, that naturally explains its strong intermittency. 

We suggest a theoretical scheme for explanation of the 
passive scalar statistics, based on smallness of the pas- 
sive scalar correlation length along the wall outside the 
diffusive layer. The scheme enables one to find explicit 
expressions for scaling exponents characterizing different 
objects. A comparison of the theoretical predictions with 
numerics shows that they agree satisfactory. Some pre- 
liminary results of our work were published in Ref. [26[ . 



The structure of the paper is as follows. In the second 
section we present our theoretical approach concerning 
the passive scalar dynamics and statistics in the periph- 
eral region and propose a scheme giving the scaling expo- 
nents. In the third section we explain our computational 
scheme, present computed moments of the passive scalar 
and integral passive scalar, give a description of the pas- 
sive scalar tongues, and compare our numerical results 
with theory. In Conclusion we give an outline of our re- 
sults, their possible applications and directions of future 
investigations. 



II. THEORETICAL DESCRIPTION 

We consider the passive scalar statistics in peripheral 
regions of a random flow, that are regions near bound- 
aries (walls). Our principal example is the viscous (lam- 
inar) boundary layer of the developed high-Reynolds hy- 
drodynamic turbulence (see, e.g., the book Q). However, 
our approach can be applied to other situations. For ex- 
ample, one can think about the peripheral region of the 
elastic turbulence [l2j. The only feature relevant for us 
is smoothness of the flow in the boundary layer, whereas 
the velocity varies randomly in time there. 

We designate the passive scalar field as 9. It can rep- 
resent both, temperature variations or concentration of 
pollutants. The passive scalar evolution (decay) in an ex- 
ternal flow is described by the advection-diffusion equa- 
tion 



d t 9 + v\70 = kV 2 6>, 



(1) 



where v is the flow velocity and k is the diffusion (ther- 
modiffusion) coefficient. Below, the fluid is assumed to be 
incompressible (that is the flow is divergentless, Vu = 0). 
The equation ([T]) implies that there are no sources of the 
passive scalar in the volume. However, we do not exclude 
a passive scalar flux from the vessel walls. 

The equation ([T]) has to be supplemented by boundary 
conditions at the wall. If 9 is the density of pollutants 
and the wall is impenetrable for the pollutants then the 
gradient of 9 in the direction perpendicular to the wall 
is zero near the wall, that corresponds to zero pollutant 
flux to the boundary. In the case we deal with the passive 
scalar decay, leading ultimately to its homogeneous dis- 
tribution in space. If 9 is temperature then its gradient in 
the direction perpendicular to the wall can be non-zero, 
that corresponds to finite heat flux through the boundary 
(from the wall) . If walls are made of a well heat conduct- 
ing material then a value of 9 (temperature) has to be 
regarded as fixed at the boundary. 

We assume that the Peclet or the Schmidt number is 
large (that is the diffusion coefficient n is small in com- 
parison with the fluid kinematic viscosity v). Then, as 



it was demonstrated in the work [22[, the passive scalar 
dynamics in the peripheral region is slow in comparison 
with the velocity dynamics. Therefore the passive scalar 
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is rapidly mixed in bulk (for a time that can be esti- 
mated as an inverse Lyapunov exponent of the flow) that 
leads to a homogeneous spacial distribution of the pas- 
sive scalar, 9 = const. The subsequent passive scalar 
evolution is related mainly to the peripheral regions that 
supply bulk by passive scalar fluctuations. We assume 
that bulk can be treated as a big reservoir, then the bulk 
homogeneous value of 9, 6b, can be regarded as indepen- 
dent of time. Below, we imply that the passive scalar 
field is shifted by 9^, that leads to the condition 9^0 
as we pass from the periphery to bulk. 

A. Correlation Functions 

Statistical properties of the passive scalar can be de- 
scribed in terms of its correlation functions 

F n {t,r u ...,r n ) = (e(t,n)...e(t,r n )), (2) 

where angular brackets mean averaging over large times 
(larger than the velocity correlation time). Since the 
velocity tends to zero at approaching the wall and the 
molecular diffusion is assumed to be weak, the passive 
scalar dynamics, determined by an interplay of the ad- 
vection and diffusion, is slower than the velocity dynam- 
ics in the peripheral region. Therefore at investigating 
the passive scalar dynamics the velocity can be regarded 
as short correlated in time and closed equations can be 
derived for the passive scalar correlation functions (see, 

e.g., EUl) 

n 

d t F n = kJ2 V ™ F " 

m— 1 

n 

+ ^2 ^2d ma [D a p(r m ,r k )d k pF n ] , (3) 

m.k—l a/3 

where the object D a p is expressed via the pair velocity 
correlation function as 

D a p(r u r 2 ) = I dt' (v a (t + t',n)v (O,r 2 )}. (4) 
Jo 

Here, again, angular brackets mean averaging over times 
larger than the velocity correlation time. 

The structure of the equation ([3]) is quite transparent: 
the evolution of the passive scalar correlation functions is 
determined by the molecular diffusion (the first term in 
the right-hand side) and by the eddy diffusion (the second 
term in the right-hand side). Therefore the quantity D a p 
can be called the eddy diffusion tensor, since it describes 
diffusion of the passive scalar related to the random flow. 
This effect can be compared with the turbulent diffusion 
of the passive scalar in turbulent flows in bulk on scales 
larger than the viscous length. However, in our case the 
eddy diffusion tensor D a p is associated with a smooth 
flow, and can be used for description of the passive scalar 
dynamics on scales smaller than the turbulent viscous 
length. 



We assume that the walls of the vessel are smooth and 
that the boundary layer width is much less than the cur- 
vature radii of the wall. Then it can be treated as flat 
in the main approximation. Let us introduce a refer- 
ence system with the Z-axis perpendicular to the wall 
and assume that the fluid occupies the region z > 0. 
Smoothness of the velocity leads to the proportionality 
laws v x , v y tx z and v z cx z 2 for the velocity components 
along and perpendicular to the wall, respectively. The 
laws are consequences of the velocity smoothness, of the 
non-slipping boundary condition v = at the wall, and 
of the incompressibility condition Vu = 0. 

Below we assume that the velocity statistics is homoge- 
nous in time, and also assume its homogeneity along the 
wall. Due to the assumed homogeneity, velocity correla- 
tion functions are dependent on time differences and on 
differences of the coordinates x and y. Say, the eddy dif- 
fusion tensor ^ is independent of time and does depend 
on differences x\ — x 2 and y\ — y 2 . However, it depends 
on both z\ and z 2 due to the strong inhomogeneity of the 
system in the direction perpendicular to the wall. A z- 
dcpcndcnce of the eddy diffusion tensor components can 
be found directly from the proportionality laws v x , v y cx z 
and v z (x z 2 . Say, 

D zz {x,y, zi\x,y, z 2 ) = fiz 2 z 2 , (5) 

where [i is a constant characterizing strength of the ve- 
locity fluctuations in the peripheral region. 

The equation for the first moment of 9 (average value 
of the passive scalar field) , (9) , is 

d t (0) =d z [fiz 4 d z (9)] + nd 2 (9), (6) 

as it follows from Eqs. ()3I5|) . Comparing the advection 
and the diffusion terms in the equation ^ one finds a 
characteristic diffusion length ru defined as 

r bl = (k/h) 1 ^. ( 7 ) 

The quantity determines the thickness of the diffusion 
boundary layer formed near the wall. Due to smallncss 
of k (remind that the Peclet or the Schmidt number is 
assumed to be large) the diffusion length is much less 
than the thickness of the peripheral region, where the 
law v z oc z 2 is satisfied. 

We consider an advanced stage of the passive scalar 
decay or a statistically stationary situation caused by the 
permanent passive scalar flux through the wall. Then 
the passive scalar 9 is non-zero primarily in the diffusive 
boundary layer, at z < ru (remind that after subtraction 
of its bulk value the field 9 should tends to zero in bulk 
that is at z — > oo). However, we are interested mainly in 
the passive scalar transport through the region z» ru, 
where the passive scalar is carrying from the diffusive 
boundary layer to bulk. There it is possible to neglect 
the molecular diffusion term in Eq. © and we arrive at 
the proportionality law 

(9) cx z~\ (8) 
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that gives the decaying rate of the average 9 as z grows. 
Note that the law (|8) corresponds to a constant average 
passive scalar flow through the planes z = const, that is 
the flux is independent of z. 

One anticipates that at z > ru high passive scalar 
moments possess some scaling behavior as well as (9). 
We introduce the corresponding scaling exponents 



oc z 



(9) 



If the molecular diffusion is irrelevant outside the dif- 
fusion boundary layer then r] n = 3 22 1. However, our 



numerical data imply that the molecular diffusion is rel- 
evant even at z > r^i (below we give an explanation of 
the phenomenon). Therefore the exponents r\ n are not 
equal to 3 and their values are subject of a special inves- 
tigation. 

Let us turn to the passive scalar correlation functions 
([2]). At z > rji their dependence on the coordinates 
along the wall are characterized by a correlation length 
I, that can be found by balance of the molecular and the 
eddy diffusion along the wall. The eddy diffusivity term 
in Eq. ((3]) can be estimated as /iz 2 , see Eq. ([5]), the 
law oc z 2 follows from the z-dependence of the velocity 
components. Comparing the molecular diffusion term 
~ k/1 2 and the eddy diffusion term in Eq. ([3]), one finds 



The quantity is of the order of r&; at z 
ishes oc z^ 1 as z grows. 



(10) 

Tu and dimin- 



B. Integral Passive Scalar 

To exclude effects of the molecular diffusion, we intro- 
duce an integral of the passive scalar field along a surface 
parallel to the wall 



Q(t, z) = A~ L J dx dy 6(t, x, y, z), (ff) 

where A is the area of the surface and z is its separation 
from the wall. We designate correlation functions of the 
object as 

*«(< ) = (Q(t,z 1 )...Q(t,z n )}. (12) 

Integrating Eq. ([3]) over x k and yk, one obtains 

<9 t $„ = re ^ + dxx ... dx n dyi ...dy n 

m=l 

d ™ \ D zzd kz F n ] + d ™ \{dkzD zz )F n ] I , (13) 

71. k— 1 m^k 

where D zz = D zz (r m ,rk). At deriving Eq. (|13l) we have 
taken some integrals in part and used the constraint 



d 

dxk 



d d 
D zx {r m ,r k ) + ——D zy {r m ,r k ) + ——D zz (r m ,r k ) = (), 

oy k oz k 



which is a consequence of the incomprcssibility condition 
= 0. 

One expects a scaling behavior of the correlation func- 
tions (H2J) at z » r bl . Then the last term in Eq. (jig)) 
can be estimated as /iz 4 9 2 $. Consequently, the term 
with the molecular diffusion k in Eq. (|13p can be ne- 
glected in the region. The argumentation is the same 
as was used for the moment (9), where the law ((SJ) was 
derived from Eq. ((6j). Note also, that in the decaying 
case the time derivative in Eq. (fTU| can be estimated as 
re/r 2 ; = y^njl, the term is much less than [iz 4 d 2 ~ /iz 2 at 
z 3> ru . Therefore the term with the time derivative can 
be neglected in the region as well, and we obtain quasi- 
stationary equations for $„ . That reflects adiabaticity of 
the passive scalar statistics. 

It is reasonably to assume that the correlation function 
F n (t, ri, . . . , r n ) is correlated along the X — Y plane on 
distances of the order of the correlation length I (|T0|) that 
is much smaller than the characteristic velocity length 
(width of the peripheral region), then D zz in Eq. (fT3j) 
can be substituted by iiz 2 n z 2 . Then we obtain closed 
equations for the correlation functions 



rn,k—l x ' 

Y, lOr ( Z m^n) , (14) 

oz m 

where we omitted the molecular diffusion term (see the 
above argumentation). 

The equations (Q3J) lead to the following closed equa- 
tions for the moments of the integral passive scalar 

dt(Q n ) = n [z 4 d 2 + 4nz 3 d z + 4n(n - l)z 2 } (6™) . (15) 

In the stationary (or quasi-stationary) case (where 
dt(Q n ) is negligible) we obtain a homogeneous differen- 
tial equation for the n-th moment that admits a power 
solution 



(9™) oc z^" 



(16) 



The exponents Cn can be easily found from Eq. (|15p . 
they are 



C„ = 2n - 1/2 + y/2n+l/4. 



(17) 



We have chosen the positive sign of the square root lead- 
ing to the reference value £i = 3, as it should be in accor- 
dance with Eq. ©. We observe anomalous scaling, that 
is a non-linear dependence of £ ra on n, that can be com- 
pared with the anomalous scaling of the passive scalar in 
the Kraichnan model [l4Hll| . 

A natural conjecture that enables one to relate the 
moments of the passive scalar 9 and those of the integral 
passive scalar O is in using the passive scalar correlation 
length / as a recalculation factor. Then we come to an 
estimation 



(6 r 



j(d-l)(n-l) 



(18) 
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Here d is dimensionality of space, that is equal to 3 in 
real flows but can be arbitrary in numerical simulations. 
The estimation ([IB]) together with Eq. (fT0|) lead to the 
relation 

Vn = Cn-(n-l)(d-l), (19) 
between the exponents introduced by Eqs. (|9I16[) . 

III. SIMULATIONS 

We conducted Lagrangian simulations where dynamics 
of a large number of particles subjected to flow advection 
and Langevin forces (producing diffusion) is examined. 
The set of the particles is used instead of the passive 
scalar field 9, that can be interpreted as density of the 
particles. A big advantage of the approach is its appli- 
cability to different space dimensions d. Indeed, in the 
scheme the number of variables (coordinates of the par- 
ticles) grows as a power of d (at a given number of the 
particles), but not exponentially. 

To establish principal qualitative features of the pas- 
sive scalar transport, we perform mainly 2d simulations. 
The setup is periodic in x (coordinate along the wall) 
with a period L. In majority of simulations the velocity 
field is chosen to be 



2irx 2nx\ L 

= z [ £i cos —j— + 4 2 sin —j— I — , 

. 2irx 2irx\ 
z I Ci sin — - 6 cos — I , 



(20) 
(21) 



where £i and £2 are independent random functions of 
time. Let us stress, that the velocity field (|2QI2ip satis- 
fies the incomprcssibility condition d x v x + d z v z = for 
any functions £1 (t) and £2 (t) ■ They are chosen to possess 
identical Gaussian probability distributions, then the ve- 
locity field (|2QI2ip has statistics homogeneous in x (along 
the wall). In our numerics, we have chosen L = 10. 

Though in reality the proportionality laws v x (x z and 
v z oc z 2 are satisfied inside the laminar boundary layer, in 
our setup the expressions for the velocity components like 
(|20I21[) are formally used at all z, that is bulk corresponds 
to z = 00. However, the law v z cx z 2 implies that the 
particles may reach the z-infinity for a finite time. That 
ensures a finite particle flux to bulk since sometimes a 
finite number of particles passages there. Therefore the 
passive scalar transport to bulk is well defined in our 
setup. 

Since the velocity correlation time in the peripheral re- 
gion is much less than the passive scalar mixing time, one 
should consider £i(i) and &(t) as white noises. However, 
in computer simulations one cannot realize zero corre- 
lation time. We model the functions by telegraph pro- 
cesses, where both functions, £1 and £2, remain constants 
inside time slots of a small duration r, and the values of 
£1 and £2 inside the slots are chosen to be random vari- 
ables with identical normal distributions. An example of 



FIG. 1: An example of the telegraph process. 



such telegraph process is plotted in Fig. [T] In our sim- 
ulations the averages were = = 1 and different 
slot sizes were used, r = 0.001; 0.002; 0.004. Then, in 
accordance with the definition (|5|), [i = r/2. 

In our scheme a particle trajectory g(t) obeys the equa- 
tion 



d t g = v(t,g) + C(t), 



(22) 



where the first term represents the particle advection and 
the second term represents the Langevin force. Let us 
stress that the variables £ are independent for different 
particles whereas the variables £1 and £2 are identical 
for all particles, according to physical meaning of the 
variables. The Langevin force C is al so modeled by a 
telegraph process with the same time slot duration r and 
with normal distributions of the values in the slots. To 
ensure a given value of the diffusion coefficient k, one 
should choose 



(£) = (C 2 Z ) = 2k/t. 



(23) 



In majority of simulations we have chosen k = r/2, and 
therefore the diffusive length was ru = 1, in accordance 
with the definition (|7J). 

Inside a time slot all the variables, £1, £2 and £, are 
time-independent constants and the equation (|22[) be- 
comes an autonomous ordinary differential equation. It 
was solved as follows. A time slot was divided into a num- 
ber of time intervals and the equation was solved (with- 
out the Langevin force) using the second order Runge- 
Kutta method. The number of intervals is z-dependent, 
it is proportional to z at z > 2.5. For z > 12 we solved 
equations for l/g z instead of g z . Both features are mo- 
tivated by the strong dependence of the velocity on z, 
v z ot z 2 . After solving the equation inside a slot a term 
produced by the Langevin force was added. To exam- 
ine role of diffusion outside the diffusive boundary layer, 
in some simulations we switched off the diffusion (the 
Langevin forces) at distances z > z^ (where a choice of 
Zd is different for different cases). 
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The particles are permanently injected near the wall in 
random positions at the beginning of each time slot. The 
simulations were performed in the stripe < z < 100, 
the particles crossing the lines z = and z = 100 were 
excluded from the set. The number of particles leaving 
the region < z < 100 through the wall is much larger 
than the number of particles escaping through the line 
z = 100. The last ones correspond to the passive scalar 
transport to bulk. A balance between the particle injec- 
tion and losses leads to a statistical equilibrium achieved 
gradually in the numerics. Thus, our simulations cover 
the statistically stationary passive scalar transport. It 
corresponds both to the steady temperature distribution 
supported by a constant heat flux from the wall and to 
the decay of the concentration of pollutants that can be 
treated adiabatically. 

An extension of our scheme to higher dimensions, d > 
2, is straightforward. We utilize the same equation (|2"2")l 
where all quantities have d components. The Langevin 
forces C are determined by the same relations (f2U| and 
a generalization of the expressions (|2QI2ip is as follows. 
The velocity v is determined by a set of 2(d— l) 2 random 
variables £nj and &ij '■ 



n—l / 
v t = z ( £uj cos 

n-l 



1 -+£ Mj Bm- F ±)-, (24) 



L 

2ixx 



L 



£ 2 jj cos ■ 



L J 7T 

2ttx 



, (25) 



where the subscripts i, j numerate first d— 1 space coor- 
dinates and the last d-th coordinate is z. Here, all ^uj 
and ^2ij are, again, telegraph processes with the same 
statistical properties as above, and we use the second- 
order Runge-Kutta scheme inside a time slot. 

In terms of the particles, the passive scalar field 8 is 
defined as a number of particles per unit volume, that is 
as a number of particles inside a box divided by the box 
volume. Of course, the definition works well provided the 
box is small (in comparison with all characteristic scales 
of the problem) and the number of particles inside the 
box is large. To satisfy these contradictory conditions one 
should deal with a large enough total number of particles. 
That is why the injection rate in our numerics is chosen 
to produce a large number of particles, 10 5 -j- 10 6 , in the 
statistical equilibrium. 



A. Tongues 

Our simulations shown that the passive scalar trans- 
port to bulk is related to specific structures of the pas- 
sive scalar. The passive scalar is concentrated mainly in 
the narrow diffusive layer near the wall. However, some- 
times a fluid jet is generated carrying the passive scalar 
from the wall towards bulk, that produces a passive scalar 
tongue with width (cross-section) diminishing as z grows. 
The property is a consequence of the law v z cx z 2 reading 




FIG. 2: An example of a passive scalar structure formed near 
the wall in the 2d random flow. Different particles are desig- 
nated by small crosses. 



that the z-component of the tongue velocity increases as 
z grows. Thus we come to a geometrical interpretation 
of the passive scalar correlation length I: it is the charac- 
teristic size of the tongue cross-section (taken along the 
wall). The cross-section behavior corresponds to an ex- 
pected decrease of the correlation length as z grows. Let 
us stress that in accordance with Eq. (TIU)) the character- 
istic tongue cross section is dependent on the diffusion 
coefficient re. 

A tongue is typically pulled from a "bump" of the pas- 
sive scalar distribution. After some time the tongue is 
tilted and then pressed back to the diffusive layer. Then 
next tongue is pulled, usually from the bump remain- 
ing at bottom of the previous tongue, and is, in turn, 
pressed back to the diffusive layer. As a result, a compli- 
cated multi-fold structure is formed, an example of such 
structure is drawn in Fig. [2J that represents a snapshot 
generated in our simulations. 

Sometimes the tongue is pulled upto z-infinity, and 
then a portion of the passive scalar (a number of parti- 
cles) is pushed to bulk. After that the tongue is tilted 
and the passive scalar current to bulk stops. That im- 
plies that the passive scalar flux, J dx 9v z in 2d, is highly 
intermittent quantity at z > r^- The conclusion is con- 
firmed by the flux histograms drawn in Fig.[3]for different 
z. In simulations, the passive scalar flux was measured 
as a number of particles crossing the plain z = const 
during a time r. At z = the flux probability distribu- 
tion is practically Gaussian, being formed by a balance 
between the random injection of the particles and their 
leaving the wall. However, the distribution becomes less 
and less Gaussian as z grows. The histograms in Fig. §3§ 
are practically symmetric. The property is related to the 
fact that only a small amount of particles in the tongue 
are pulled to bulk, the majority of particles returns back, 
that produces practically equal fluxes to bulk and to- 
ward the wall. That explains why at z > ru the root 
mean-square fluctuation of the flux is much larger than 
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FIG. 3: Histograms of the passive scalar flux at different sepa- 
rations from the wall. The root mean square fluctuations are 
much larger than the average value and the histograms are 
practically symmetric. At z = the probability distribution 
is Gaussian whereas at z > ru it has exponential tails. 
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FIG. 5: Log-log plot of the moments of 9g, (0%), multiplied by 
z 3 , for 8 = 0.03125 and n — 14-6 in the case where diffusion is 
substituted by a constant velocity carrying the particles from 
the wall. 
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FIG. 4: Log- log plot of the moments of 6s, (&s), multiplied 
by z 3 , for S = 0.03125 and n = 1 4 6. The graph reflects sim- 
ulations where diffusion occurs everywhere, and is switched 
off at z — 3 or z = 12. 



its average value. 



B. Moments 

Based on numerical data one can compute moments 
and correlation functions of different quantities charac- 
terizing the passive scalar statistics. One can consider 
both local functions or integral objects. All the quanti- 
ties are computed as time averages. 

We introduce an object 6s that is a number of particles 
inside a square box of size 5 divided by its area S 2 (in 
2d). The quantity 9g is close to 9 provided the number of 
particles is large and the size of the box is small enough. 
Moments M n of the quantity 6s, M n = for n = 

1 4- 6 are computed as averages over time intervals 10 6 4- 



10 7 r as S = 0.03125. The results are presented in Fig. [4] 
where the moments multiplied by z 3 are plotted in log-log 
coordinates (solid curves). We see that the prediction flSJ) 
for the first moment is perfectly satisfied whereas higher 
moments deviate strongly from the diffusionless law oc 
z~ 3 . We conclude that the diffusion is relevant at z > ru, 
indeed. 

To verify the conclusion we repeated the simulations 
switching off diffusion at z > 3 and at z > 12. Results 
are presented in Fig. 2] by dashed curves. We see an ap- 
pearance of plateaus, starting just from z = 3 or z = 12 
and corresponding to the law cx z~ 3 , in accordance with 
Rcf. [22| . The plateaus are observed in restricted re- 
gions of separations from the wall z slightly diminishing 
as n grows. An explanation of the fact is that cutoffs 
of the plateaus are observed where 5 becomes of the or- 
der of the passive scalar correlation length (along the 
wall). To check this conjecture we repeated the simula- 
tions for larger values of 5 and observed that the plateaus 
shrink as S grows. That confirms our explanation. To 
be absolutely sure that the diffusion is relevant we per- 
formed simulations without diffusion but with adding a 
constant velocity V carrying the particles away the wall. 
Results of the simulations are presented in Fig. [5] where 
the moments (Of), multiplied by z 3 , are plotted. We see 
plateaus signalling that outside the boundary layer the 
passive scalar moments behave in accordance with the 
diffusionless prediction. 

The next object of our investigation is the integral 
quantity that is the passive scalar integrated along 
the wall, see the definition (fTTj) . In numerics it is de- 
termined by a number of particles in a slice of thickness 
5, parallel to the wall, divided by its volume (area), we 
designate the ratio as <ds ■ In our 2d setup the area is 
equal to LS, where 5 is chosen to be much less than z. 
The moments of Q$, (&g), are computed by time aver- 
aging over a long time ~ 10 7 r. To check robustness of 
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FIG. 6: Moments of in log-log coordinates, n = 1 -¥■ 6. In 
the region z > ru the results collapse onto single curves for 
three times r = 0.001, 0.002, 0.004 and four different values of 
the diffusion coefficient. 
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FIG. 7: Moments of <d$ in log- log coordinates, n = 1-^6. The 
results are obtained for two cases where the diffusion occurs 
everywhere and where it is switched off at z > 3, and also for 
two different velocity fields: with two and four harmonics. 



the results we performed computations for different time 
slots, t = 0.001, 0.002, 0.004, and for four different values 
of the diffusion coefficient k. The figure |6] demonstrates 
that the values of each moment collapse to a single curve 
in the logarithmic coordinates \n(z/ru) and \n((<d n ) / C n ) 
where the factors C n are the corresponding moments near 
the wall. 

One can check that, in accordance with our theoret- 
ical expectations, the moments of Qg are insensitive to 
diffusion. To illustrate the assertion we present in Fig. [7] 
the moments of Qg computed at r = 0.002 for two cases: 
in the first case the diffusion occurs everywhere and in 
the second case it is switched off at z > 3. One can ob- 
serve no difference between the data. One can also try 
to use an alternative to the expression (|2QI2ip velocity 
field. There are two series of data plotted in Fig. [7] and 
corresponding to the velocity (|20I21[) and to the velocity 
with four random factors 



TTV X 2ttx 2ttx 

—— = £i cos — — + ii sin — — 
Lz L L 

4nx Attx 
+£ 3 cos — + £4 sin — , 

v z 2ttx 2ttx 
-2=6 Sln S2 cos — 

AlTX AltX 

+2£ 3 sin — 2£ 4 cos — — . 



(26) 



(27) 



Again, there is no visible difference between the series. 

We observe that the moments of O are decreasing func- 
tions of z that are power- like in the region z > ru- Ex- 
tracting the scaling exponents Cn, sec the definition (|17|) . 
for n = 1 -7- 6 in 2d we obtain values that are presented in 
Fig. |S] as the lower set of points (some smooth curve is 
drawn through the points for better visualization). We 
conducted analogous simulations for higher dimensions, 
upto d = 5. The results are depicted in the same figure 
[SI We see that the exponents Cn depend on d, however, 
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FIG. 8: Exponents of the moments (Og), for n = 1 -r 6 and 
space dimensions d = 2 -j- 5. For comparison the theoretical 
curve Cn = 2n — 1/2 + \/2n + 1/4 is plotted (solid line). 



for d > 3 they are close to the theoretical values (fTT)) 
represented by a solid line. 

One can think that the deviations from the theoreti- 
cal values (fTT)) are related to an existence of additional 
passive scalar (relatively long) correlations along the wall 
that can be produced by the multi-fold structures of the 
type drawn in Fig. f2J The long correlations should lead 
to increasing moments of the passive scalar in comparison 
with the short correlated case. It is naturally to expect 
that the fold effect becomes less pronounced in higher di- 
mensions. Indeed, Fig. |S| shows that the deviations from 
the values (|17p diminish as the space dimensionality d 
grows. That confirms our explanation. 

To check our conjecture, we conduct simulations for 
the velocity field containing, like the expressions (|26I2T[) , 
a set of harmonics in terms of the period L: the ninth, 
tenth and eleventh ones. In such velocity field corre- 
lations related to the multi-fold tongue structures have 
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FIG. 9: Exponents of the moments (GJ), for n = 1 -f- 6 
and space dimensions d = 2 for two different velocity fields: 
containing only first harmonic and three (ninth, tenth and 
eleventh) harmonics. For comparison the theoretical curve 
Cn = 2n - 1/2 + y/2n + 1/4 is plotted (solid line). 
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FIG. 10: The difference of the scaling exponents of moments 
for the integral passive scalar and for the passive scalar, C„ — 
?7n, computed at 8 = 0.03125 in 2d. For comparison, the 
theoretical prediction n — 1 is drawn. 



to be suppressed, and, consequently, the exponents £ n 
should be close to the theoretical values (jl7[) . This ex- 
pectation is confirmed by our simulations, the results are 
presented in Fig. [SJ where measured exponents arc plot- 
ted. Wc sec a good agreement of the measured and the- 
oretical exponents. 

One can extract from our numerical data the exponents 
i] n of the moments of 9s, see the definition (|9|), as well 
as Cn ■ It is interesting to check the theoretical prediction 
(fTTJ)) . For the purpose we plotted the difference — i] n 
as a function of n, see Fig. [TO] For comparison, the 
theoretical straight line n — 1 (correct in 2d) is drawn in 
the same figure. We see a good agreement, confirming 
the scaling (fTU|) of the passive scalar correlation length 
l(z) along the wall. 



IV. CONCLUSION 



We performed extensive numerical simulations of the 
passive scalar mixing in peripheral regions of random 
Sows like high-Reynolds turbulence. The simulations 
confirm earlier theoretical expectations and reveal a lot of 
new details. At advanced stages of the passive scalar mix- 
ing passive scalar fluctuations are concentrated mainly in 
the narrow diffusive layer near the boundary provided the 
Pcclct or the Schmidt number is large. We found that 
the passive scalar transport from the diffusive boundary 
layer to bulk is related to passive scalar tongues formed 
by jets directed to bulk. The tongues arc objects respon- 
sible for strong intermittency characteristic of the passive 
scalar transport through the peripheral region. 

We examined the passive scalar statistics outside the 
diffusive boundary layer and realized that both, moments 
of the passive scalar 9 and of the passive scalar inte- 
grated along the wall, O, possess well pronounced scal- 
ing in terms of the separation from the wall z. We com- 
pared the corresponding exponents extracted from our 
numerics with our theoretical scheme and established 
their agreement. However, one should be careful since 
our theoretical predictions are correct for an infinite ves- 
sel and can be violated for the numerics where the ve- 
locity correlation length along the wall coincide with the 
velocity period. We found also an agreement between the 
theoretical prediction for the tongue cross-section depen- 
dence on z and our numerics. Therefore the simulations 
confirm our theoretical predictions. 

There remain some problems to be solved in future. 
We are going to extend our consideration incorporating 
average flows (like in pipes) which arc shear-like near the 
wall. Another natural extension of our approach is re- 
lated to chemical reactions in random flows. One can also 
note polymer solutions, where the polymer elongation is 
very sensitive to the character of the flow. The problem 
is significant, e.g., for the elastic turbulence. However, a 
long-time memory characteristic of the polymer solutions 
could modify our results. We considered smooth walls in 
our work. There is a set of questions related to rough- 
ness of the wall, possible corners, caverns and peaks. All 
the objects can modify our conclusions, it is a subject of 
special investigation. 

Our results agree qualitatively with data known from 
investigations of turbulent plumes in turbulent flows 
where a complicated space structure of the passive scalar 
fluctuations is observed [27T - l3lj ]. We believe that statis- 
tical properties of the structure can be explained on the 
basis of our results implying production of the passive 
scalar tongues pushing to bulk. The explanation needs 
a generalization of our scheme where turbulent velocity 
fluctuations in bulk should be included. That is a subject 
of future work. 
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